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We study the effect of strong magnetic fields, of the order of 10 18 — 10 19 G, on the instability 
region of nuclear matter at subsaturation densities. Relativistic nuclear models both with constant 
couplings and with density dependent parameters are considered. It is shown that a strong magnetic 
field can have large effects on the instability regions giving rise to bands of instability and wider 
unstable regions. As a consequence we predict larger transition densities at the inner edge of the 
crust of compact stars with strong magnetic field. The direction of instability gives rise to a very 
strong distillation effect if the last Landau level is only partially filled. However, for almost completed 
Landau levels an anti-distillation effect may occur. 

PACS numbers: 21.65.-f 26.60.Kp 26.60.-c 97.60.Jd 



I. INTRODUCTION 

Recent investigations seem to show that soft 7-ray repeaters and some anomalous X-ray pulsars are neutron stars 
which may have surface magnetic fields larger that 10 15 G [1—3], the so called magnetars. Until recently, the strongest 
estimated magnetic field is of the order of B = 2 x 10 15 G and was detected in a quite young star, SGR 1806-20 [4]. 
According to [5] a fraction as high as 10% of the neutron star population could be magnetars. 

The effect of the strong magnetic fields on the equation of state (EOS) of stellar matter in neutron stars has been 
studied both at low densities below neutron drip, of interest for the study of the outer crust of neutron stars [6] , and 
at high densities, of interest for the study of the interior of compact stars [7, 8]. In this last case field-theoretical 
descriptions based on the non-linear Walecka model (NLWM) [9] were used and several parametrizations compared. 
It was shown that they have an overall similar behavior. It was recently shown in [10] that the EOS at subsaturation 
densities, including densities of the order of the densities at the inner edge of the crust of a compact star, was 
particularly affected by fields of the order of 10 18 G. 

An important characteristic of the nuclear matter is the appearance of a liquid-gas phase transition at subsaturation 
densities. The role of the isospin is of particular importance. Indeed, since nuclear matter is composed of two different 
fluids, namely protons and neutrons, the liquid-gas phase transition can lead to an isospin distillation phenomenon [11]. 
The region of instability is determined by the spinodal curve. Due to the symmetry energy, the EOS of /^-equilibrium of 
magnetic free nuclear matter is thermodynamically stable. The stability of the EOS is determined from the curvature 
of the free-energy: a positive curvature corresponds to thermodynamically stable matter. 

If we consider stellar matter at very low densities, nuclei in matter are expected to form a Coulomb lattice embedded 
in the neutron-electron sea that minimizes the Coulomb interaction energy. With an increase of the density, nuclear 
"pasta" structures emerge [12]. The existence of pasta phases may modify some important processes by changing 
the hydrodynamic properties and the neutrino opacity in supernova matter and in the matter of newly born neutron 
stars [13]. Also, the pasta phases may influence neutron star quakes and pulsar glitches via the change of mechanical 
properties of the crust matter [14]. It is therefore important to study how the magnetic field could affect the extension 
of the pasta phase and the isospin distillation effect. 

In fact, sufficiently strong magnetic field affect the extension of the unstable region. In order to have a better 
understanding of the effect of the magnetic field on the instabilities of nuclear matter at subsaturation densities we 
study in the present work the effect of a strong magnetic field on the thermodynamical spinodal instabilities obtained 
from the free energy curvature matrix [15, 16]. Recently, it was shown that the magnetic field and Joule heating have 
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the important effect of maintaining compact stars warm for a longer time [17]. This kind of simulations need the EOS 
of the crust. It is, therefore, important to make a study that shows when should the magnetic field be taken explicitly 
into account in the EOS of the crust. An unstable region in a wider density range will correspond to a larger crust 
and the properties of the star depending on the crust will be affected. 

In the present paper, we will consider two relativistic effective approaches: a NLWM, TM1 [18], with constant 
coupling parameters, and a density dependent relativistic hadronic (DDRH) model TW [19] with density-dependent 
coupling parameters. DDRH models seem to give results closer Skyrme interactions than NLWM, at subsaturation 
densities [20] . 

In section II we make a brief review of the models, the EOS under the effect of a magnetic field and the stability 
conditions. Results arc discussed in section III and conclusions are drawn in the last section. 



II. THE FORMALISM 



A. EOS of nuclear matter under a strong magnetic field 

In the present section we make a short review of the field-theoretical approach used to obtain the EOS of nuclear 
matter. Within this approach, the baryons interact via the exchange of a, uj and p mesons in the presence of a uniform 
magnetic field B along the z-axis. We start from the Lagrangian density of TW [19, 21] model 

£= C b + C m . (1) 

b— n,p 

The baryon (b=n, p) and meson (a, lo and p) Lagrangians are given by 



A = *b [il^ - gi7„A" -m b + r a <T - I^w" - ^T p T^^ - ^p N K b a^F^ * b , 
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- -^F^ + tmlp^-tp^P^ (2) 



respectively, where are the baryon Dirac fields. The nucleon mass and isospin projection for the protons and 
neutrons are denoted by nib and r 3 b = ±1, respectively. The mesonic and electromagnetic field strength tensors are 
given by their usual expressions: f2 M „ = dpLu u — P pv = d p p v — d^p^, and F^ v — d p A v — d v A p . The nucleon 

anomalous magnetic moments are introduced via the coupling of the baryons to the electromagnetic field tensor with 
G\xv = \ b7t)7i>] and strength Kb with n n = g n /2 = —1.91315 for the neutron and k p — (g p /2 — 1) = 1.79285 for the 
proton, respectively, and where gi are the Lande g factors for the particle i (g p — 5.5856912 and g n = —3.8260837), 
and /i at is the nuclear magneton. The electromagnetic field is assumed to be externally generated (and thus has 
no associated field equation), and only frozen- field configurations will be considered. The electromagnetic couplings 
are denoted by q. The parameters of the model are the nucleon mass nib — 939 MeV, the masses of mesons ni a , 
m u and m p and the density dependent coupling parameters T which are adjusted in order to reproduce some of 
the nuclear matter bulk properties and DBHF (Dirac Brueckner Hartree-Fock) calculations [22], using the following 
paramctrization 



where x — p/p sa t, with 



and, 



Ti(p) = Ti(p sat )fi(x), i = <r,u,p (3) 



1 + bj (x + djf 

fi{x)=<H- ; —2, i = (T,u, (4) 

1 + Ct (x + di) 



f p (x) = exp [-a p (x - 1)] , (5) 



with the values of the parameters raj, Tj, a^, c, and di, i — cr,oj,p given in Table I for TW model [19]. Other 
possibilities for these parameters are also found in the literature [23] . 
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TABLE I: Parameters of the TW model. 



The symmetry energy and its first and second derivatives are important to understand the instability region. NLWM 
models become very stiff above saturation densities while DDRH models have a softer behavior. On the other hand, at 
subsaturation densities DDRH models have larger symmetry energies and a larger extension of the instability region 
for very asymmetric matter. In Fig. 1 we compare the symmetry energy of the models we will consider in the present 
study: TM1 and TW. As expected TM1 has a smaller symmetry energy at subsaturation densities and a larger one 
above the saturation density. 




P b (fm 3 ) 

FIG. 1: Symmetry energy for all models used in the present work. 

Notice that in the DDRHM the nonlinear meson terms are not present, in contrast with the usual NLWM. For 
TM1 model we add to the Langrangian density, Eq. (21) , with gi = Tj, 

C = \bm n {g„af + \c{g a af + if S J( S i/) J . 



The coupling parameters are constant and given in [18]. 

From now we take the standard mean-field theory (MFT) approach and display only some of the equations. A 
complete set of equations and description of the method can be found in the literature {e.g., Ref. [7, 24]). For the 
description of the system, we need the baryonic density, the energy density of nuclear matter, and the pressure. The 
energy density of nuclear matter is given by 

e= J2 £b + \ m ^ a ' 2 + \ m2 ^ + \ m lpl ( 6 ) 

b— p,n 
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where the energy densities of nucleons have the following forms 
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For the neutrons we have introduced 

m n = to* — sfiNKnB, 

where the effective baryon masses are given by 

to* = to — T a <7. 
The pressure of the system is obtained from the expression 

Pm = ^2 Pbpb - £■ 
b—n,p 

The energy spectra for protons are neutrons are given by 

El 



yjm* 2 + 2vq p B — s/j,nK p B 

( 1 1 

- 7; s ^NH n B + -m n 



(7) 

(8) 
(9) 

(10) 

(11) 
(12) 

respectively, where v = n + | — sign(q)^ = 0, 1, 2, . . . enumerates the Landau levels of the fermions with electric 
charge q, the quantum number s is +1 for spin up and —1 for spin down cases, and the rearrangement term is given 

by 

&£, 
dp 
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where p s — p P + p s n and pt, = p p + p n , with the expressions of the scalar and vector densities for protons and neutrons 
given by [7] 
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where k v F v s and k F s are the Fermi momenta of protons and neutrons which are related to the Fermi energies E P F and 
E F through 
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,p2 F p2 

K F,u,s — F 



kf s = Ef-ml 



^m* 2 + 2vq p B — sp^K p B 



The chemical potentials of nucleons within TW are given by 



,i p = E P F + IW + ^I>° + S* 



(15) 



(16) 



For TM1 we have similar expressions with the last term, the rearrangement term, equal to zero. 
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B. Stability conditions for nuclear matter 

At subsaturation densities nuclear matter has a liquid-gas phase transition and homogeneous matter is not stable 
within a given range of densities. The stability conditions for asymmetric nuclear matter, keeping volume and 
temperature constant, are obtained from the free energy density T, imposing that this function is a convex function 
of the densities p p and p n . For stable homogeneous matter, the symmetric matrix with the elements [25, 26], 



T%i — 



d 2 T 

dptdpj 



(17) 



known as stability matrix, must be positive. This corresponds to imposing that the trace and the determinant of Tij 
are positive. In terms of the proton and neutron chemical potentials pi, the stability matrix is given by 



(18) 
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with m = f£| 



OJ- I 

dpi \ T >PJ^- 

For nuclear matter, the largest eigenvalue of the stability matrix is always positive and the other becomes negative 
at subsaturation densities. We define the thcrmodynamical spinodal at T=0 as the curve on the p n , p p plane defined 
by the points for which the determinant of T%j is zero; that is, the smallest eigenvalue is zero. Inside the region limited 
by the thermodynamical spinodal the smallest eigenvalue of J-ij is negative and nuclear matter is unstable [26] . At 
T = 0, T is equal to the energy density defined by Eq. (6). The eigenvalues of the stability matrix are given by 



A± = \ (Tr(^) ± v/Tr(^)2 - 4Dct(J0) . 



(19) 



The stability condition requires that they are both positive. When one curvature becomes negative the system is 
thermodynamically unstable and can decrease its free energy by going in the instability direction [26] , defined by the 
direction of the eigenvector associated to the negative eigenvalue. The eigenvectors 5p± of the stability matrix are 
given by 



\± — j-j 



i,j =p,n. 



(20) 



In the following we study the direction of instability inside the spinodal section for both models considered. 

C. Stability conditions for npe matter 

Stellar matter at low densities is formed by protons, neutrons and electrons in equilibrium with respect to weak 
interaction processes. Until now we have considered the subsaturation instability region of neutron-proton (np) nuclear 
matter. In this section we investigate the effect of the inclusion of electrons on the stability conditions of nuclear 
matter when electrons are included. In particular, we will calculate the thermodynamical spinodal sections for npe 
(neutron-proton-electron) neutral matter. Since matter is neutral the proton and electron densities must be equal, 
i.e. p p = p e . 

Electrons are included in a minimal way in the system, and are described by the following the Lagrangian density 



C e = V> e {ilpd* 1 - q e lpA^ - m e ) ipe 



(21) 



where ipi are the electron Dirac fields and m e = 0.511 McV. 
Including the electrons, the stability matrix (18) becomes 
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(22) 
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and the stability conditions arc equivalent to the ones indicated in the previous subsection: the trace and the deter- 
minant of T must be positive. 
The electron density is given by 



\qe\B 
2tt 2 



(23) 



where hp v s is the electron Fermi momentum related to the Fermi energy Ep by 

kf^ s = Ef-(ml + 2v\q e \B). 

For npe neutral matter 

T = e + e e 

where e was defined in (6) and the electron contribution is given by 



_ \0e\B x ^ 



k l p^ s Ep + (ml + 2v\q e \B) In 



^ml + 2u\q e \B 



(24) 



(25) 



We next discuss np nuclear matter and npe neutral matter, /^-equilibrium stellar matter is a particular case of npe 
neutral matter, with the proton and electron fractions defined by chemical equilibrium conditions, namely 



with Up and a n defined in (16) and u e = Ep. 



III. RESULTS AND DISCUSSION 



In the present section we first show the dependence of the spinodal section on the magnetic field, both for neutron- 
proton {np) and neutron-proton-electron {npe) matter. From the crossing of the /3-equilibrium EOS of stellar matter 
with the thermodynamical spinodal we make a prevision of the transition density and transition pressure at the inner 
edge of the crust of a compact star. 

For a zero magnetic field, the direction of instability of nuclear matter gives rise to a distillation effect, corresponding 
to the formation of droplets of dense matter with low isospin asymmetry in a background of a neutron gas with a 
small fraction of protons. This effect has been observed experimentally in heavy-ion reactions [11]. Therefore, we also 
discuss the effect of the magnetic field on the direction of instability, namely in which way it affects the distillation 
effect. 



A. Spinodal section np matter 



We will first consider np nuclear matter and determine the instability region limited by the spinodal surface. In 
Fig. 2 and Fig. 3 we show the spinodal sections obtained making A_ = 0, where A_ was defined in Eq. (19), on the 
{p p , p n ) plane for TM1 and TW and for several values of the magnetic field. We define the magnetic field in units of 
the critical field 3% — 4.414 x 10 13 G, so that B = B* B%. For a field with the intensity B% the electron cyclotron 
energy is equal to the electron mass. 

We present the numerical results both not including and including the contribution of the anomalous magnetic 
moment (AMM). In all figures where both cases are considered we show on the left panel the results without the 
magnetic field and on the right panel the results including the AMM. 

The magnetic field has a strong effect not only on the size, giving rise to larger instability regions, but also on the 
shape of the instability zones which is not symmetric with respect to the p p = p n line, contrary to np matter for 
B = 0. The increase of the instability region is due to Landau quantization which softens the EOS. For magnetic 
fields B* > 2 x 10 5 , the protons are totally polarized for all the values of the densities considered and the size of the 
spinodal zone is larger than the one obtained for B = 0. Including AMM decreases the spinodal region with respect 
to the results without AMM for all the values of the magnetic field considered. This behavior is explained by the 
extra stiffness that the inclusion of AMM brings into the system due to neutron spin polarization. 
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FIG. 2: (Color online) Spinodal section on the p p , p n plane for TM1 at T = MeV and for severals values of magnetic fields, 
B — B* B e , a) without and b) with AMM. For B* = 10 5 the spinodal section is formed by two separate regions. 
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FIG. 3: (Color online) Spinodal section in terms of p p versus p n for TW at T = MeV and for severals values of magnetic 
fields a) without and b) with AMM. For B* = 10 5 and 2 x 10 5 the spinodal section has respectively three and two separated 
parts. 



The effect of the Landau quantization on the spinodal section is explicitly seen in the spinodal for B* = 10 5 for 
TM1 [dashed line Fig. 2 a)]. The spinodal region consists of two separate zones each one corresponding to one Landau 
level, the one corresponding to the first Landau level extends to larger neutron densities than the one corresponding 
to the second level. In order to understand this effect, we plot in Fig. 4 the proton chemical potential for B* = 0, 
B* = 10 5 and B* = 3 x 10 5 . We have identified the instable regions with thick lines. It is seen that for the larger 
field the proton chemical potential changes smoothly because for all the densities shown only the first Landau level 
(LL) is occupied. At low densities the chemical potential decreases with density and only above 0.025 fm -3 it starts 
increasing with density. For B* — 10 5 the proton chemical potential shows a cusp corresponding to the end of the 
first LL and beginning of the second LL. The unstable regions, in this case, correspond to the beginning of each 
LL when the slope of the chemical potential is smaller. The smaller the magnetic field the larger the number of LL 
occupied at subsaturation densities and the larger of independent sections which make up the whole spinodal section. 
For reference we include the chemical potential at B — 0: it increases smoothly with density with a quite constant 
slope. 

If the AMM is included the instabilities regions are smaller, as discussed above. The structure (bump) appearing at 
p n ~ 0.05 fm~ 3 is due to the neutron polarization: below that value of the density the neutrons are totally polarized. 
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For the TW model with B* — 2 x 10 5 and 10 5 , Fig. 3, we also get a spinodal region formed by several bands, 
respectively two and three bands. An interesting feature of this model, is that the band with the largest Landau 
level may extend to larger neutron densities than lower levels. This does not occur for NLWM and has to do with 
the behavior of the symmetry energy which increases in a smoother way for DDRHM than for NLWM above p = 0.1 
fm -3 . As a result the slope of the chemical potentials is not so large. 
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FIG. 4: (Color online) The proton chemical potential for B = and for two magnetic field intensities, B* = 10 5 and B* — 2x 10 s , 
obtained within the TMl model for the neutron density p n — 0.05 fm -3 and excluding the AMM. The thick lines represent the 
regions of instability. 

Another feature of the spinodal regions with magnetic field and without AMM, that is present in both models we 
have studied, is the extension of the spinodal for zero proton fraction: in the absence of the magnetic field there is 
no instability but the inclusion of the magnetic field changes this behavior: the instability region at p p = extends 
until a finite p n value, model dependent, but independent of B. This value is ~ 0.15 fm -3 for TW and ~ 0.186 fm -3 
for TMl. 

To understand this behavior seen at zero proton fraction, we consider the TMl model. For p p = 0, we obtain the 

fc n3 

corresponding finite value of p n = r^y from the Fermi neutron momenta kp solution of the equation T)et(J-) = 
with p p = 0. Explicitly, the latter equation can be written as follows 
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This equation is independent of the magnetic field and, therefore, all the spinodal regions for different magnetic fields, 
without AMM, in Figs. 2 and 3 have the same value of p n for p p = 0. The inclusion of AMM changes this feature: 
p n is still finite for p p = but its value is B dependent. 

For models with density dependent couplings the spinodal region extends to larger densities for the larger proton 
densities when compared with NLWM. This is due to the behavior of the symmetry energy: while for NLWM the 
symmetry energy increases quite steeply for densities above saturation densities, DDRH models have a much smoother 
behavior and the symmetry energy of these models take much smaller values than NLWM for densities above p — 0.15 
fm~ 3 . 




FIG. 5: (Color online) Proton and neutron chemical potentials as function of the proton density for p n = 0.02 fm -3 and for 
B* — 10 a) without and b) with AMM for TM1 (dashed line) and TW (full line). The thick segments of each curve represent 
the regions of instability. 



In Fig. 5 we plot p p and p n for B* = 10 5 with the models TM1 (dashed lines) and TW (full lines). The thick 
segments of the curves lines identify the instability regions defined by the two (TM1) and three (TW) bands which 
form the spinodal. These curves are obtained for a fixed neutron density, p n = 0.02 fm~ 3 . It is seen that above 
p p — 0.1 fm -3 the proton chemical potential within TW is much softer and this seems to be the reason for the 
appearance of the third band in this model. 

For each model TM1 and TW we identified the crossing density, and corresponding pressure, of the EOS of j3- 
cquilibrium stellar matter and the corresponding spinodal for each value of the magnetic field considered. The EOS 
of state was obtained considering neutrons, protons and electrons in /3-equilibrium. In order to illustrate what was 
done we represent in Fig. 6 the spinodal sections obtained within TM1 for B* = 0, 10 5 and 5 x 10 5 respectively by 
full, dashed and dotted thick lines. We include in the same figure, using thin lines with the same type of curve for 
each B value, the corresponding EOS of /3-equilibrium stellar matter. The crossing spinodal-EOS is identified by a 
big dot. Both the spinodal and the EOS are plotted in the p p , p n plane. 

The crossing density of the EOS with the thermodynamical spinodal gives a prevision of the transition density 
[27, 28] to an homogeneous phase, and is always larger than the one obtained from the crossing of the EOS with the 
dynamical spinodal for npe matter, which includes the Coulomb interaction. In [29] the authors have shown how the 
transition density and respective pressure were related to the fraction of the star's moment of inertia contained in the 
solid crust, and obtained a relation between the radius and mass of compact stars. 

In Tables II and III the values of the crossing density and respective pressure are given for stellar matter under 
different magnetic fields, respectively without and with the AMM. The values of the crossing density for B = 0, 0.069 
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FIG. 6: (Color online) Spinodal section in the p p , p n plane for TM1 at T — MeV and for severals values of the magnetic 
field without AMM. For each value of the magnetic field, it also plotted the EOS of stellar matter in /3-equilibrium (thin lines). 
The crossing of the EOS with the respective spinodal, large dot in each spinodal, represent the transition density to continuous 
matter. 



TABLE II: Predicted density, proton fraction and pressure at the inner edge of the crust of a compact star at zero tempera- 
ture, as defined by the crossing between the thermodynamical instability region of np matter and the /3-equilibrium EOS for 
homogeneous, neutrino-free stellar matter in the p P ,pn plane. The AMM is not included. 
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0.30219 


1.3795 


5 x 10 s 


TM1 


0.16783 


0.40921 


1.6324 




TW 


0.19784 


0.40194 


2.3310 



fm 3 for TM1 and 0.085 fm 3 for TW, can be compared with the corresponding ones obtained from the crossing 
of the dynamical spinodal with the EOS Ref. [20, 27], respectively 0.06 fm" 3 for TM1 and 0.075 fm- 3 for TW. As 
expected they are a bit larger, with TW model having a larger crossing density than the other. The effect of the 
magnetic field is to increase the values of the crossing density: at B* = 10 5 both models have similar transition 
densities of the order of ~ 0.1 fm -3 corresponding to a much larger pressure for TM1 than TW. For B* = 3 x 10 5 
the transition densities increase to p ~ 0.14 — 0.15 fm~ 3 . 

In Table III we show the same data given in Table II but including the AMM in the calculation. The conclusions 
are similar: the transition density increases with the increase of the magnitude of the magnetic field but not so fast. 
However, the corresponding pressures are larger than before. We conclude that the existence of a strong magnetic 
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TABLE III: Predicted density, proton fraction and pressure at the inner edge of the crust of a compact star at zero temperature, 
as defined by the crossing between the thermodynamical instability region of np matter and the /3-equilibrium condition for 
homogeneous, neutrino-free stellar matter. The case where the AMM is included. 



Models 



Y„ 



P m { MeV fm" 



10 5 

2 x 10 5 

3 x 10 5 
5 x 10 5 



TM1 
TW 
TM1 
TW 
TM1 
TW 
TM1 
TW 



0.086942 
0.091391 
0.096337 
0.092438 
0.11046 
0.11251 
0.12822 
0.14880 



0.17670 
0.17829 
0.29468 
0.30512 
0.37135 
0.38035 
0.47093 
0.48803 



1.3801 
1.2809 
1.5373 
1.3549 
1.7091 
1.8047 
1.6616 
2.7717 



field at the crust gives rise to a larger crust. 



B. Spinodal section npe matter 



We have studied the effect of the magnetic field on the instability region of np matter in the previous section. For 
npe matter without magnetic field, NLWM models still present a small thermodynamical instability region but for 
DDRHM models there is no instability region [15]. The incompressibility of the free electron gas is so high that the 
spinodal disappears or almost disappears. 

In Fig. 7 the spinodals for npe matter are shown for TM1 and different magnetic fields. In fact, although including 
electrons, the instability region can become almost as large as the B = np-spinodal. This is due both to the Landau 
quantization of the orbital motion of protons and electrons: the incompressibility of the electron gas is smaller than 
the one of a magnetic free electron gas. 
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FIG. 7: (Color online) Spinodal section in terms of p p versus p n for TM1 for npe neutral matter at T = MeV and for severals 
values of magnetic fields (a) without and (b) with AMM. 



For TW, contrary to the B — case, the inclusion of the magnetic field gives rise to a spinodal region as seen in 
Fig. 8. The behavior of this model with the magnetic field is similar to TM1. We also point out that the inclusion 
of the AMM has a strong effect on the spinodal part corresponding to the first LL: it is drastically reduced or even 
disappears. 
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FIG. 8: (Color online) Spinodal section in terms of pp versus pn for TW for npe neutral matter at T = MeV and for severals 
values of magnetic fields (a) without and (b) with AMM. 



C. Direction of instability 



The eigenvector associated with the negative eigenvalue of the free energy curvature matrix defines the direction of 
the instability and tells us how does the system separate into a dense liquid and a gas phase. It was shown in [15, 20] 
that in the absence of the magnetic field the direction of instability favors the reduction of the isospin asymmetry 
of the dense clusters of the system, and increases the isospin asymmetry of the gas surrounding the clusters, the so 
called distillation effect. This effect is represented in Fig. 9 where it is seen that for the B ~ curve (thick full line) 
the fraction Sp~ /Sp~ is larger than p p j p n below y p = 0.5 and the other way round above. 

In Fig. 9 we plot, for TM1, the ratio 5p~ /Sp~ for p = 0.06 fm -3 as a function of the proton fraction. Several 
results for different values of the magnetic field are shown by the thick lines. The thin lines represent the ratio p p j ' p n 
for reference and y p = 0.5 points corresponding to symmetric matter, as well as Sp~ /5p~ = 1, which is the ratio of 
density fluctuations for symmetric matter with no field. For the largest field considered the spinodal region contains 
a single Landau level and the curve varies smoothly starting at 8p~ /5p~ ~ 1.5. We point out the very large value of 
this fraction, always above 1, for y p < 0.5. The magnetic field favors an increase of the proton fraction quite above the 
symmetric matter value. For B* = 10 5 the spinodal has two bands, see Fig. 2 and 3 , corresponding to the occupation 
of the first two Landau levels. The transition from one to the other is clearly seen with a large discontinuity of 
5p~ /5p~ at y p ~ 0.7. Above this y p value the curve behaves like the previous ones. However for y p < 0.7 the behavior 
is quite different: the curve decreases from the value at y p =0, which is independent of the magnitude of the magnetic 
field, to a value much smaller than the corresponding value of the fraction p p j p n . The fluctuations will not drive the 
system out of the first Landau level and therefore the larger the proton fraction, the closer the system comes to the 
top of the band and the smaller are the allowed proton fluctuations. For y p > 0.7 or for the larger magnetic fields 
the Landau levels are only partially filled and the fluctuations will never drive the system out of the corresponding 
Landau level. 

Similar features are obtained for TW and/or including the AMM. In Fig. 10 we show, respectively for TM1 (top) 
and TW (bottom), the fraction Sp~ /Sp~ as a function of y p for a fixed baryonic density, p = 0.06 fm~ 3 , chosen inside 
the instability region. For y p > 0.5, TM1 and TW behave in a similar way, while below this value the main difference 
is the smaller 6p~ /Sp~ for TW, corresponding to a smaller distillation effect. This behavior is also present for B = 
and it was shown that this was due to the presence of the rearrangement term. The inclusion of the AMM favors 
larger proton fractions because neutron polarization stiffens the EOS. 

In Figs. 11 and 12 we represent the fraction 5p~ /Sp~ as a function of density respectively for two values of y p , 0.2 
and 0.4, for np matter (thick lines) and npe matter (thin lines). We consider TM1 and TW. Both models have a very 
similar behavior for finite values of B although for y p = 0.2 and B = they differ: for TM1, 5p~ /Sp~ increases with 
density while, for TW, this fraction decreases for p > 0.02 fm~ 3 . This effect is not so strong for y p = 0.4 and for npe 
matter the fraction is always quite small due to the presence of electrons which prevents large proton variations. 



13 




For B — 10 5 , Sp~ /5p~ decreases with density while for B = 5 x 10 5 the opposite occurs. In both cases only the first 
Landau level is occupied, however for the lower field the first Landau level is almost full and the density fluctuations 
will occur in such a way that the system stays in the same Landau level: the larger the total density the smaller 
the fluctuations. For the larger field the first Landau level is only partially filled, far the top of the band. For the 
same nuclear density, the larger the proton fraction the lower the system energy and therefore the fraction Sp~ /Sp~ 
increases with density. 

In Fig. 12 we give the same information with y p = 0.4. While for B* = 5 x 10 5 for the range of densities considered, 
matter occupies only one Landau level, for B* = 10 5 we may have two (TM1) or three (TW) Landau levels, see Figs. 2 
and 3. This explains the discontinuities occurring for ~ 0.1 fm -3 . For the highest magnetic field only one Landau 
level partially filled comes into play and therefore the fraction Sp~ /5p~ increases with density because that is favored 
energetically. For B* = 10 5 the presence of almost filled Landau levels prevents the existence of large proton density 
fluctuations. 



IV. CONCLUSIONS AND OUTLOOK 



In the present work we have studied the instabilities of np matter and npe neutral matter under very strong 
magnetic fields. The fields considered are much stronger than the strongest field measured until now at the surface 
of a magnetar which is B* ~ 10 2 for SGR 1806-20 [4]. However, it is expected that fields in the interior of neutrons 
stars will be much larger. The present work shows how fields of the order of B = 5 x 10 18 G could affect the inner 
crust of a compact star. 

We have considered two relativistic nuclear models: one NLW model (TM1) and one DDRH model (TW). For both 
models, we have determined the spinodal surface, from the curvature matrix of the free energy, for different magnitudes 
of the magnetic field. It was shown that the instability region could be divided into several bands according to the 
magnitude of the magnetic field and the number of the Landau levels occupied. The presence of the magnetic field 
will generally increase the instability region. Making a crude estimation of the transition density at the inner crust 
of a compact star under a strong magnetic field from the crossing of the EOS of /3-equilibrium stellar matter with 
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FIG. 10: (Color online) &p~ /6p„ plotted as a function of the proton fraction with T = MeV and p = 0.06 fm -3 for the TM1 
(top) and TW (bottom) and for severals values of the magnetic fields without (a) and (c) and with (b) and (d) AMM. The 
fraction p P /p n is given by the thin dotted line. 



the thermodynamical spinodal, it was shown that the transition density and associated pressure increases with the 
magnetic field. This will affect the structure of the star increasing the fraction of mass and of the star' s moment 
of inertia concentrated at the crust. These effects will be noticeable if, for densities of the order of 0.1 fm -3 , the 
magnetic field is of the order of B* = 10 4 or larger. 

The TW model has larger instability regions than the TM1 model for the larger proton densities. A smoother 
increase of the proton chemical potential for the first model justifies this behavior. This behavior of the symmetry 
energy may even give rise to a larger number of bands in the spinodal of TW than the spinodal of TM1 for the same 
magnetic field. 

We have also investigated the direction of instability. It was shown that if the Landau level is only partially occupied 
the density fluctuations are such that the system evolves for a state with dense clusters very proton rich immersed in 
a proton poor gas. A larger proton fraction is favored energetically due to the degeneracy of the Landau levels. If on 
the other hand, we study the fluctuation of particles occupying an almost complete Landau level, proton fluctuations 
cannot be so large and it may even occur an anti-distillation effect with a decrease of the proton fraction of the dense 
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FIG. 11: (Color online) 8p p /5p n plotted as a function of the density with T — MeV and y p — 0.2 with (thin lines) and 
without (thick lines) electrons for the TMl and TW models and for severals values of the magnetic fields (a) and (c) without 
and (b) and (d) with AMM. 



clusters. This is due to the fact that these fluctuations will keep the system in the same Landau level. 
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